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Abstract 

The effect of random velocity field on the kinetics of single-species and two-species annihi- 
lation reactions is analysed near two dimensions in the framework of the field-theoretic renor- 
malisation group. Fluctuations of particle density are modeled within the approach of Doi. The 
random incompressible velocity field is generated by stochastically forced Navier-Stokes equation 
in which thermal fluctuations - relevant below two dimensions - are taken into account. 

1 Introduction 

The effect of density fluctuations on the asymptotics of reaction rates in low dimensionalities has 
attracted considerable attention recently ||l], |2|. Reaction rates may also be affected by fluctuations 
of an advective velocity field. Most work in this direction has been carried out for the case of 
quenched random drift [0, ||, |^. Recently, the asymptotic behaviour of the unimolecular reaction 
A + A —> (!) in a dynamically generated random drift has been analysed with the aid of field-theoretic 
renormalisation group (RG) 0. In my report I shall describe a similar approach to the bimolecular 
reaction A + i? — > 0. 

Classical rate equations for densities ua, ns with the homogeneous initial condition n^(0) = uao, 
ns(0) = HBO 

—— = -KqUaUb , —J- = -KqUaUb (1) 
dt dt 

yield the "normal" decay laws. Let, for definiteness, uao < ^so, then 

nAO ^-KoinBQ-nAo)t 



UAit) ^ (nBo-^Ao)— e-^«("^™)*, nsit) ^ {ubo - uao) ■ (2) 



In case of equal initial densities uao = nBo = no a powerlike decay takes over: 
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n^(t) = nB{t) ^ — — . (3) 



Kot 



A heuristic account of initial-density fluctuations amplified by diffusion together with numerical 
simulations ||^, || yield anomalous decay at < 4, e.g. for uao = ^bq = Uq: 

which is slower than the "classical" decay rate @. 



2 Second quantization for reactions 

For a systematic analysis of the effect of density fluctuations it is convenient to use a field-theoretic 
approach. There are two possibilities available for this. First, the more widely known Martin-Siggia- 
Rose (MSR) approach |||], which corresponds to the solution of a Langevin-type stochastic equation 
to describe fiuctuations. Physically, this approach is well-suited to situations in which fiuctuations 
are of "external" origin, e.g. small-scale thermal fiuctuations from the point of view of macroscopic 
(hydrodynamic) scale physics, or fiuctuations caused by some external random source. Second, the 
several times reinvented Doi approach |[T^, in which the randomness is described by a probability 
distribution functional (PDF) or, if you like, an infinite set of probability distribution functions on 
a lattice, with the subsequent set of master equations instead of Langevin equations. This approach 
is better suited to cases, in which fiuctuations of intrinsic origin are dealt with. This is the case, for 
instance, for density fiuctuations due to randomness in the (chemical) reaction process itself. 

Therefore, I will use here the Doi approach, which allows for a "mesoscopic" analysis of density 
fiuctuations. To calculate expectation values with the probability distribution functional (PDF) 
P({n^(x)}; {nB(x)},t) for the particle densities n^(x) and ^^(x), the formal solution of the set of 
master equations for the PDF may be expressed in a functional form with the aid of bosonic field 
operators p!0| , pT| with the commutation relations 



[V^A(x),V^i(x')] =5(x-xO, [M^),M^')] = [^l(x),^+(x')] =0 



and a similar set for the B particles. In a fairly standard-looking second-quantization setting the 
average of an observable O may be written as a vacuum expactation value 

m)) = E 0[{n^},W}]P(W)};W},t) 

{ni(x)} 

= {0\O[{^+ + l)^A,{^+ + l)^B]e-^'' (5) 

with the kinetic operator 

H' = Jd^ |^lV(v^^) + ^t^iv^B) - Dao^I^^^a - DBoi^i^^i^B 

+ Ko{jP^ + ijt + i^A^i)^AM. (6) 
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The last exponential in (|5|) corresponds to the initial PDF. A customary choice is the Poisson dis- 
tribution for the local particle number ||l|, However, in view of the hostile nature of interaction 



I have allowed for negative initial correlations by using in (|^ a bivariate Poisson distribution ||T2 
This choice leads to the following expressions for the low-order moments of initial densities: 



nA(x, 0) = uao , ?^b(x, 0) = hbo , AnA(x, 0)AraB(x', 0) = -ro^/nMn^S^,^', 

AnA(x, 0)AnA(x', 0) = uao 5x,x' , Anij(x, 0)Anij(x', 0) = ubo 5x,x'- (7) 

Physically this corresponds to thermal fluctuations with anticorrelations in initial reactant densities. 

3 Dynamic action for the advection-diffusion-controlled reac- 
tion A + B ^ 

Construction of perturbation theory through the T exponent for the evolution operator 

allows to write the expectation value (|5|) as the following functional integral: 

(0(t)) = J V[i,] 0^(1, 1, ^3^'^'' , (8) 
where Ojv is the normal symbol of the operator O 

5*1 is the dynamic action ||^ 

- DAoi^lv^i^A - Dsoi^t^^i^B + Ko {^pt + + i^ti^s) i^Ai^B } (9) 

and 5*0 contains terms brought about by the initial bivariate Poisson distribution 

So = Jdx Hao^a + ^Boipt - roV^AoUBQipA^t ■ 
Schwinger equations with respect to ipAi i^B 



dt^A + V(v^a) - DaoV^^a + Ko^a^b) = nAoS{t) , (10) 
dti^B + V(v^b) - DBo^^i^B + Koi^Ai^B) = riBoSit) , (11) 



in the mean-field approximation yield the classical rate equations (|T|) for the homegeneous average 
densities (nyi(t)) = {ipAit)), (^^(t)) = (ipBit))- It should be borne in mind, however, that the second 
and higher order moments of the fields ipA and i/jb are not equal to the corresponding moments of 
the densities. For instance, the pair correlations of A particles are given by (n^(t, x)nyi(t, x')) = 

{[ijA{t, ^)Mt, X') + 5(X - ^')Mt, X)]). 
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To describe fluctuations of the drift fleld v in (H) I use random velocity fleld generated by the 
transverse stochastic Navier-Stokes equation 

dtv + P(v ■ V)v - uoV\ = r (12) 

with the incompressibihty conditions: V ■ v = 0, V ■ f = 0. For the random force the Gaussian 
distribution with zero mean and the correlation function 

(/^(xi,t0/:(x2,t2)) = 5(ti-t2) / ^P.„(k)dKA;)e^'-(^-^^) (13) 

is assumed. In (|T3D Pm„(k) = 5mn — kmkn/k'^ is the transverse projection operator in the wave- vector 
space, and df{k) is a function of the wave number k and the parameters of energy pumping, which 
is used to produce stationary random drift. The kernel function is often chosen in the nonlocal form 

df{k) = gi.ulk^-''-^^ (14) 



to generate turbulent velocity field with Kolmogorov's scaling ||T3|, (which is achieved by choosing 

6 = 2). 

The stochastic Navier-Stokes equation (0) is a Langevin type equation leading to a MSR operator 
functional which in the functional-integral form gives rise to the following action functional 



+ / (itdx V 



5*2 = - y^tfixfix' v(x, t) ■ v(x', t)c?/(|x — x'l) 

-dt^ - (v • V)v + z/qV^v] . (15) 

Combined averaging over density and velocity fluctuations yields 

{0{t)) = J V[^, v] 0^(1, ^A, 1, i^B) e^«+^^+^^ (16) 

for the expectation value of the observable O. 



4 Decay asymptotics controlled by stable fixed points 

Power counting shows that in the case in which all three reaction terms in (||) have the same scaling 
dimension the critical dimension of the model is two Q. Near two dimensions, however, the drift 
part ([T5|) of the dynamic action with the nonlocal kernel ( |T4|) is not multiplicatively renormalizable. 
Therefore I have used the kernel function (l^ with a local term added at the outset: 

dfik) = g.^vlk^-"'^^ + g^o^lk^ • (17) 

Apart from rendering the fleld theory multiplicatively renormalizable, the local term also has an 
important physical meaning: with a suitable choice of the parameter it describes thermal fluctu- 
ations of the velocity fleld near equilibrium. 
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Taking this into account, I write the renormalised action in the form 



S 



- UAyZ2Ai^A^'^i^A - UByZ2Bi^B^'^i^B 



— -V 

2 

+ V- 



(18) 



+ dx 



9tv + (v- V)v-z/ZiVV]} 

KA/i'^2'5A^A + K^B^-^Z^Bi^B " PAi'^^sV'a V^B 



in which, apart from the standard renormalisation of the dynamic action, also the renormalisation 
of the initial condition - predicted by power counting and confirmed by calculations - is introduced. 

Renormalisation constants have been calculated in one-loop approximation with the use of com- 
bined dimensional and analytic regularisation with the parameters e and 5 = {d — 2)/2, which 
eventually give rise to a two-parameter expansion of critical exponents and other physical quantities. 

The unrenormalized parameters of the initial conditions have positive canonical scaling dimen- 
sions. This means that the corresponding running parameters grow in the long-time large-scale 
limit. Therefore, some kind of partial summation of the perturbative expansion is called for to cope 
with this problem. A natural way would be the use of skeleton equations for Green functions with 
dressed field averages and correlation functions. In the case of single-species annihilation reaction 
A + A ^ ^ this leads to well-controlled estimates of the behaviour of scaling functions in the long- 
time limit [|l|. Basically, this amounts to independence of the asymptotics of the initial density ||16| . 
In the present case of bimolecular annihilation, however, a similar direct summation has not been 
found 0. Unfortunately, I have not been able to do any better with an analytic solution of the set 
of integro-differential equations, which can be written for the dressed one-point and two-point Green 
functions of the present model and incorporate the effect of initial conditions completely. 

In the leading order in the coupling constant A the initial density fluctuations change the classical 
rate eqs (|5) by the addition of an inhomogeneous term. This leads to the system 



dtUA = 

dtriB = 



. -25 i Pf^'^ 1 
-Az//i <^ uaUb - J- — J ■ -Tjj^ \ , 

X -25 \ PP'^ \ 

-^""^ ["'""'-[AMua + ubWA^ 



(19) 



with the initial condition: nyi(O) = kaP'^Z^a, ^b(O) = ^bP'^Z^b- For the important special case 
ua = ub = n the special Riccati's equation results. The solution is known and may be expressed 
in terms of modified Bessel functions -ft'2/(4-d) and ^/(A-d)- In this solution the blow-up of initial 
conditions is controllable and the asymptotic behaviour of the particle density may be inferred. 
However, the effect of initial conditions is not fully accounted for by the manageable system of eqs 
( |T^ ) , and therefore the following results give the correct asymptotic density decay with this provision 
only. 

The renormalised action (^) gives rise to a system of characteristic equations with four IR stable 
physical fixed points with the following asymptotic decay of the density. 
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(i) Gaussian fixed point 

gl = g; = X* = 0. 

The Gaussian fixed point is stable, when e < 0, 5 > 0. Asymptotic decay of the number density in 
terms of physical (unrenormalised) parameters 

is not given by the classical solution (|) but is slower. This is different from the unimolecular case @ 
in which at the Gaussian fixed point the mean-field solution holds. Note that there is no dependence 
on the rate coefficient in (pO|), but the parameter of initial correlations remains. The infiuence of 
initial correlations becomes irrelevant and the mean-field decay oc 1/t is restored only at > 4 |||]. 

(ii) Thermal fixed point 



gl = 0, g; = -327c6, u* = ^^^^ ^ , A* = -27r(v^ - 1)5 . 

The basin of attraction of this fixed point is 5 < 0, 2e + 35 < 0. Decay rate is faster than the 
initial-density-fiuctuation induced: 

„M . f ;V'\ (21) 



Here, r is a reference time scale and 5 = (i/2 — 1. 
Reactive kinetic fixed point 



64 7re(2e + 35) , 647r 



9x = —r — , 92 



9 e + 5 9 5 + e 



17 — 1 An , 

is stable, when e > 0, — |e < 6 < — |e. To linear order in 6, e, decay exponent the same as in thermal 
fixed point 

^ ,iT\..„. r,y'- (22) 

f— >oo 



The independence of the exponent of time of 6 and e in (|2lD and (|22D is most probably an artifact 
of the one-loop approximation. 
(iv) Passive kinetic fixed point 



, 64 TT 6(26 + 35) , 647r , Vl7 - 1 

3i = ^ — ' ^2 = ^T— , u = , A =0, 

9 6+0 9o+6 2 

is stable, when 6 > 0, 5 > — 16. Decay rate is faster than the initial-density-fiuctuation induced here, 
too: 

t—>oo 



[47rz/o(VT7 - l)r]'^/4 Vt 



Here, r is the reference time scale. The decay exponent in (^) is exact. 
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Figure 1: Basins of attraction of the Gaussian fixed point (G), thermal fixed point (T), the reactive 
kinetic fixed (Kl) and the passive kinetic fixed point (K2) in the (5, e) plane. 



From these results it follows that the decay exponent is a continuous function oi 5 = d/2 — 1 and 
e - apart from logarithmic corrections on the basin boundaries. 

As in the case of unimolecular reaction P|, the fixed point corresponding to the pure diffusion- 
limited reaction 

gl = gl = Q, A* = -47r(v^-l)5 

is unstable m d < 2. This means that any velocity fluctuations (including the ubiquitous thermal 
noise!) drive the system to the advection-diffusion-controlled regime with different decay exponents. 

In case of unequal initial densities the system ([T^ ) leads to a particular form of Riccati's general 
equation, for which the solution seems to be unknown. Therefore, I will not discuss this case here 
in the hope that a solution of this equation in terms of known special functions may eventually be 
found with a reasonable effort. 



5 Conclusions 

The main physical result of this work is that diffusion-limited two-species annihilation reaction is 
shown to be unstable to short-range (thermal) velocity fluctuations for d < dc = 2 and unstable to 
long-range (turbulent) velocity fluctuations for d > 2. Decay exponents in four stable advection- 
diffusion-controlled regimes have been calculated at one-loop order. Although the renormalisation 
and fixed-point analysis are fairly straightforward, the asymptotic analysis of scaling functions is 
not. The blow-up of the effective (running) initial conditions in the scaling functions is difficult 
to control, in contrast with the single-species case, and a firm conclusion about the asymptotics of 
scaling functions is not yet available. 

As to possible generalisations, it would be interesting to amend the decay analysis by a similar 
treatment of a stationary state with reactant sources. 
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It seems quite plausible that the asymptotic behaviour of the density heavily depends on the 
localisation of the initial density profiles; therefore, an analysis of the problem with localised initial 
conditions would be desirable. 

Due to the incompressibility condition imposed on the drift field, the present results have a direct 
physical meaning at d > 2 only. Thus, the effect of compressibility should be analysed. This, 
however, does not seem to be feasible at present in the full stochastic Navier-Stokes framework. 
Therefore, to make some progress in this direction, it would be interesting to analyse the effect of 
velocity fiuctuations with given statistics instead of dynamically generated random drift. 
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